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PREFACE 


Matrix algebra has been used extensively in analysis 
of passive electrical circuits and grid-driven vacuum tube - 
amplifiers, and thus has forged our intuition for the worr- 
ings of passive reciprocal devices and active unilateral 
devices. The advent of the transistor, which is not entire- 
ly unilateral | has caused the need of wider knowledge of 
Non-unilateral active circuits, which include feedback amp- 
lifiers. Examples of the above types of circuits are shown 
in Figure 1. 

This paper presents methods of matrix algebra applied 
to the stability problem. The analysis must lead to the 
Bame result as any other method, and, although the amount 
of calculation may be the same as for other methods, the 
advantages of matrix algebra lie in arrangement and organ- 
ization of the problem. 

The writer wishes to thank Professor R, Kahal for 


suggesting the problen. 
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Figure 1 - Active and Passive Circuits 
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(Listed in the order of their use in the text) 


Numbers in modified parenthesis, / 7, refer to the Biplio- 
graphy. In general, matrices will be represented by pagren- 
theses, except where the omission of parentheses causes no 
ambiguity -~ in such cases capital letters represent the 
matrix, 1.e., 
2 Z 
7 11 12 

421 422 
det F is the determinant of the matrix F 
Iz| represents the absolute value of Z 


ny, is the transpose of the matrix n, 80 that the elements 


By of the matrix n equal the elements Bay Of ne. 


the amplifier transmission characteristic 
transmission characteristic through the feedback network 


the return difference 


a 8) 
| 


the return ratio 
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CHAPTER I 
» INTRODUCTION 


1. Early literature — 

The analysis of linear four-poles by the use of matrix 
algebra has received iicdvive Gudsiarnt ih the literature 
f2,5,117. The treatment, at first restricted to passive 
networks obeying the reciprocity relation, was later ex- 
tended to networks containgng vacuum tubes /3,4,15/ and gy- 
rators /197. Early literature regarding stability of net- 
works failed to make use of matrix algebra /14/. Examples 
of the application of the transmission matrix to problems 
of stability are in papers by Jefferson /O 7, Sulzer /177/, 
and Hinton /7_7, but these are limited merely to the deri- 
vation of the matrix network function corresponding to the 
Nyquist 4/7 criterion. Honnell / 8 / extended these, methods 
and wrote @ generalized criterion for stability of a closed 
tandem connection of two terminal-pair networks. More re- 
cently, Mason /137 indicated means of transforming ampli- 
fiers into stable forms and gave a relation between power 


gain and stability in feedback amplifiers. 





CHAPTER II 
STABILITY 


1. General 

An amplifier may be Said to be stable if its response 
ultimately decays to zero after it has been subjected to 
& smal) impressed disturbance which itself dies out. If 
the response increases until limited by non-linearity, the 
amplirier is said to be unstable. 

The determination of the stability of a non-unilateral 
amplifier leads directly to the analysis of the linear 
four-pole shown in Figure 2, which may include internal 


feedback loops. 





Figure 2 - The general non-unilateral amplifier 


The characteristic equations may be written 


Ei * Z1111 + Zj0lI0 2 Eg = Zely (a 
at 
Eo 2 2931, + Zoalo e Zz Jo 





or, more concisely, in matrix fora, 


& 21. 412) /11 
| 1 = | . (2) 
Eo Zo, 422/ \ta 
In general, we may say that Z)5 and Zo; will not be 
equal, since the feedback amplifier includes active elements 
which are not bilateral, although there may be included pas- 


sive, bilateral elements and active unilateral elements (1.e., 


vacuum tubes). Then, from equation (1), 


Eg(Z2e a Z) 
fi, 


7s 
(3) 
ina ~Z2 Ee 
oe oA 


We can see that, with zero impressed voltage, I; and 
I, can be other than zero when A =- 0, with A a function 
of j4. When A is considered a function of the more general 
frequency variable, p= + jW, so that the currents are of 
the form e(7 + Jj )t, it can be seen that for A = 0, stabil- 
ity of the system requires 7 <0. 

Thus, the location of the zeros of A are limited to the 
left half-plane, so that all the free modes of response are 


positively damped and so vanish with time. The presence of 





@ zero in the right half-plane is indicative of a negative- 


ly damped or runaway response. 





CHAPTER III 
VARIOUS STABILITY CRITERIA 


lL. Return-Difference and the Nyquist Diagram 

In a single loop structure the fundamental quantity 
is the loop Ererees Peto F » Which can be described as 
the voltage returned around the circuit when the loop is 
broken just behind any grid and a unit am.f. is applied 
to that grid. The factor 1 “APs which measures the re- 
duction of the effect of parameter variations in a feedback 
amplifier, is closely related to the loop trangmission. 

Bode {/ 1_/ used these quantities for a general defin- 
ition of feedback, and defined the quantities Returm Dif- 
ference and Return Ratio. The two quantities are so close- 
ly related that the term feedback may be applied to both. 

The return difference, F, is defined as the potential 
difference appearing between the grid node N and the grid G 
if this connection is broken and, while all other genera- 
tors in the network are dead, a generator of unit e.m.f. 
is inserted between cathode K and the grid G (see Figure 
3). 

The return ratio, T, is the voltage between the grid 
node N and the c,thode K and is equal to (P-1), It will 
be remembefed that this concept is applicable when opening 
the input to one grid. Thus, in an ordinary amplifier, 


4 PSP: 








Figure 3 - Definition of Return Difference 


The return difference, F, is shown to be equal to A/ Aw 
where A is the circuit determinant, and A°%1s the cir- 
cuit determinant with the grid input under question opened. 

The requirements on A for stability of the structure 
were mentioned previously as being that the roots could 
not appear in the right half-plane. Evidently we wish to 
investigate the behavior of F in the right half-plane to 
determine the stability or inatability of a circuit, but 
the actual operation of the circuit is of interest only 
at real frequencies, i.e., on the jw <-gxis. 

The Nyquist diagram is useful in an ynvestigation of 
the behavior of the circuit 4 in the right half-plane in 
that it utilizes a theorem from the theory of functions of 
@ complex variable, which states that: 


If a function f(z) 1s analytic, except for pos- 
Bible poles, within and on a given contour the 


6 





number of times the plot of f(z) encircles the 
origin of the f(z) plane in the positiye direc- 
tion, while z itself moves around the prescribed 
contour once in a clockwise direction, is equal 
to the number of poles of f(z) lying within the 
contour diminished by the number of zeros of 
f(z) within the contour, when each zero and 

pole is counted in accordance with its multi- 


plicity, / 

The concepts of return difference and return ratio 
are useful in determining the stability of a network in 
that the Nyquist diagram is essentialjy a plot of the 
values of F, a function of p, corresponding o values of 
p lying on the boundary of the right half of the p-plane. 

Practical reasoning assumes F to be negligible for 
large values of p, ee at high frequencies the circuit 
becomes capacitive (this is the case while p moves around 
the infinite semy-circle). The complete diagram then 
becomes a plot of only the real frequency values of F 
for -jJw < pc+ jw. (See Figure 4). 

For the case of a single return loop, 4° has no zeros 
in the right half-plane, and any zeros of 4 1n the right 
half-plane will result in the plot of F encircling the 
origin; thus, the condition for stability is that the plot 
shalj] not encircle the origin. 

Although the above ae soription envisions a plot of F, 
in practise the diagram is plotted i, terme of the return 
ratio, T, 60 that the critical point shifts one unit to 
the left to (-1,0) in the T diagram, or to (+1,0) in the 


Me diagram. ‘ 





ju, -. > y semicircle 
.. of infinite 
radius 





APT? 





Figure 4 - The Nyquist diagram 





2. The Return Difference Matrix 

Bode's concept of the return difference was applic- 
able to opening one grid input. This was extended to open=- 
ing several grid inputs simultaneously, by the use of matrices 
by Tasny-Tschiassny /187. He defined the return difference 
matrix: 


The return difference matrix is defined by the 
statement that the voltage matrix 

V - FE 
appears between the grid nodes N and the grids 
Gif the voltages, or currents, or ocoth, gener- 
ated by the input generators are made zero, and 
the considered elements are disconnected from 
the grid nodes N and connected to the cathodes 
K through zero impedance generators producing 
the voltage matrix E, /18 


The V's and E's are column matrices and F is a square 


matrix, so that from the definition above, 


aT ha tra ** * Ss Ey 
Vo Foy oe ee e@ @ e« @ E, 
oe a) (4) 
ue Fal © 6 e @ e@ “Fon En 


The return ratio matrix, T, determines the correspond- 
ing voltages appearing between the grid nodes N and cath- 
odes K, so that F and T are square matrices, the number of 


rows and colymns being the number of considered elements. 





Then, extending Bode'’s definition of T, 


fo Pio’ * ‘Tin Z1. Ti * Ty) [2.0 Oro 
Poy ee eee @ ; T53° e ee e i O 1 Qeeeere 
@ e = © s OQ Leercee (5) 
Fal - =” ¢ @ hae Thi’ eee a O iii 
or, more concisely, 
ye T tI (6) 


where I is the appropriate unit matrix. 

The elements of the return difference matrix are ob- 
tained by selecting, successively, voltage matrices E in 
which one element is equal to unity and all others are 
zero. The analysis must be carried out as many times as 
there are considered elements. 

The condition for the network with considered grid 
inputs opened to be stable is shown to De 

det FZ O, (7) 
where det F is the determinant of the matrix F, so that 
g > exists. | 

Finally, the criterion for stability with considered 

elements ditowed is shown to be, for det F= 0, 7<0O, 


where F is a function of the complex frequency p-=¢7 + j. 


LO 





In terms of the Nyquist diagram, points on the jw -axis 
in the p-plane are mapped into the det F plane, and the re- 
quirement for stability is that the origin not be encircled. 


3. Honnell's generalized criterion. 

Another approach to the stability criterion is shown 
by Honnell / 8/7, who wrote a generalized criterion for stab- 
Llity of a closed tandem connection of two terminal.pair 
networks. 


The transmission matrix of one two terminal.pair is 


E a & E 
( 1) ‘( 11 2),( : (8) 
Ty Qo, 859/ \nIo 
and the overall transmission matrix for n two terminal. 


pair networks connected in tandem as shown in Figure 5 is 


E E 
( ) (al).(a@). . can). - (9) 
I, «I 


n 


Or, in more concise form 


E E 
( a a) | . : (10) 
Ty -I, 
where (a°) is the overall transmission matrix, and is ob- 


tained from the matrix product of the associated at matrices 


for the individual two terminal-pairs, taken 4, the proper 


1l 





order. 





Figure 5 - Tandem connection of two terminalepairs 


The Overall stability criterion, as a function of the 


complex frequency, p, is shown to be 
A(p) = a°,+ ang ~ (det a®) - 120 (11) 


when the tandem connection is closed into a single loop. 
The roots of A(p), where p lsthe complex frequency 7+ ju, 
must lie in the left half-plane, so that the requirement 
for stability is 7< 0. 

The methods shown can establish that a circuit is 
stable or unstable, but give no indication of means of in- 
creasing the stability of a circuit on the verge of instab- 
ility. | 
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CHAPTER IV 
THE UNILATERAL AMPLIFIER 


1. General 

It is obvious that an unstable amplifier can be stabil- 
ized by utilizing any expedient that will maintain the zeros 
of O& in the left half-plane, or, in effect, alter the course 
of the Nyquist diagram and so avoid the enclosing of the 
critical point. This may be accomplished by altering the 
scale of the feedback, or by transforming the network into 
&@ form which will be absoltuely stable. With only a few 
Simple restrictions, an amplifier which is unilateral (i.e., 
Z)> Of the impedance matrix is equal to zero) 1s absolutely 
stable. 

This chapter shows that a unilateral amplifier may be 
obtained with the aid of matrix algebra, and the conditions 
under which this amplifier is stable. The actual transforn- 
ation 1s made with a coupling circuit which may be represented 
by matrices in various ways, and when the transformation has 
been made with a lossless coupling circuit, the new amplifier 
which appears at the terminals of the coupling circuit has the 
same power gain as the former. Once the amplifier is made 
unilateral, additional coupling may be considered for feed- 
back and designed according to elementary feedback theory. 

. Figure 6 illustrates how the coupling networks are con- 
nected to the amplifier. In many cases, the feedback net- 


work may be combined with the unitlateralizing network. 


Isee Appendix II for application of matrices to feedback. 
13 





The amplifier is considered to have a common ground at each 
end, and is shown as a three terminal network rather than 
868 4n Figure 2, The characteristic equations are the same 


asp for a two terminal-pair, 





ACTIVE DEVICE 







COUPLING 
NETWORK 


SOURCE 


Figure 6 - Connection of coupling networks 


2. Positive-Definiteness 
If the matrix of the amplifier considered as a two 
terminal-pair is written 


Roy + JXq] Rio + IJXy10 
Roy + JXo} Roa + JX 
it may be shown /12_/7 that the system is stable for any 


passive termination at all frequencies if the following 


14 





conditions are satisfied: 


Roo > 0 


2 
4(Ry Roo t XyoX%0,) (Ry yRoo - RyoRoq) - (RyoXaq - RoyXz0)" 7 0 


These are the Gewertz conditions for stability in which 
Peek ingedance matrix is described as being positive-definite. 
If the amplifier is unilateral, 1.e., 2,5 = 0, these 
conditions reduce to 
Ri, > 0 (19) 
Ron > O 
Thus, the conditions for stability of a unilateral amplifier 
reduce to rather simple forms, and it is for this reason 
that a means is sought to transform an amplifier which may 
be suspiciously near instability to a form which is absolute- 
ly stable. 
The following section shows that a coupling network 
can always be found which will make the transformed matrix, 
Z', unilateral. The requirements for stability, then, are 
that Rj}, and Roo are both positive. When either Ri or Roo 
is negative, a new lossless coupling can be found which 


Will make both Rj, and R455 positive. 


05 





3- Iosslees reciprocal transformations 
If the amplifier of Figure 6 is non-unilateral, it 


is characterized by the equations written at terminals a: 


Typ ta.> 
SOURCE & sb Pe 
¢ , 
| 
2 2 
| 
Ips> la. 


Flgure 7 - The amplifier and coupling 


* 


Fao “oa “pe/ \ed 
or, in more abbreviated forn, 
Bo - 4ao1, (14) 
In the same manner, the equations of the amplifier 


and coupling circuit may be measured at terminals b, 
Ey, = Zpplp (15) 
In order to transform oan into Zpp? the equations 


of the coupling network C alone are written. The relations 


in the multi-terminal network C between FE), Eo, Ij, and Ip 


Wirth othe 


16 





may be expressed by 


es = K a (16) 
Ty Io 


where the E's and I's are column matrices:such that 


Ey = (Poa) etc., 


Fal 


and Kis a 4 x 4 square matrix whose elements are independ- 


ent of IT and E, K, then, is given by the submatrices 


° 


and equation (16) is equivalent to 


Ey AE, + BI 


(17) 


Ty CE, + DI5 
The elements of A and D are dimensionless numbers, B is. 
an impedance m,trix, and C is an admittance matrix. 

The characteristic equations of the coupling network 


Can then be written in the forp 


E 
(as (3 4 
Ep Ip 


Where Z 1s an impedance matrix derived from K (167 such 
" 





that 


| i Ac-lp - 4 
z= 
c=l cwlp 


If the coupling cirguit is reactive, equation (18) 
may take the form 


Ea\. /-JCan SJCay Ig | 
(3). (Hm) a 
Eb “J%a—q = SPpp ly 
where the C's are real; or, in expanded form, 
Eg = -JCoat, + JCoplp 


En = aloe, 5 9 ICpolp 


(Equation (19) has been expanded into equation (20), but equa- 


(20) 


tion (20) still represents abbreviated matrix notation 


from equations (15), (14), and (15)). 
The transformation from Za, into Zpp May be made by 


combining equations (14), (15), and (20), so that 


(Zoo Wi JCaglle = JCaptp (21) 


(Zyp - SCyplIp =-JCpala 
from which, 
(Zpyp - ICpp) = CyalZaq + JCgg)~+ Cap (22) 


The coupling circuit C is reciprocal since we labeled 
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it reactive, and matrix C is symmetric. Thus, Cag and Chp 
may be arbitrary, but symmetric,* and C,, = (Chades where 
the subscript t represents taking the transpose. Such a 
network, then, transforms the impedance matrix measured 

at the amplifier terminals into one at the coupling tern- 


inals by operations of the general form: 
z' = jx' 2 n(Z + 4x)7* ny (23) 


where x and x’ are arbitrary real symmetric matrices, yn 
1s an arbitrary real matrix, and ny is the transpose of n. 
It can be shown that lossless reciprocal coupling can 


accomplish any of the following operations: 


Z'2oZ + jx (with x, = x) 
es Zt (24) 


Z' = nZnt 


where x and n are real. 


4, Unil@teral Gain 
Mason / 187 shows that the parameter U, the unilateral 


gain, is given by: 


2 
454 ~ 410 
. (25) 


A(R Roo ~ RyoRoy) 





and 18 invariant under the three operations of transforn- 


See Appendix I. 
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ation given in (24)for lossless reciprocal coupling. For 


example, under the real transformation 

Z* = nZnt (26) 
consider the real part of (26): 

R' = nRn, (27) 


or 


t e 
sq “2, iY fas i mal a (28) 
 ]  ] 
Roi Roe M51 "o0/ \Roy Roa} \8o1 M20 
whose determinant is 
2 
R11R40 - RyaRo, = (nyyMoo - MygMaq)“(RyyRoo ~ RypRoy) (29) 
Also transpose equation (26) and subtract, yielding 
Zz’ - 2 sn(Z - 2,)n, (30) 
whose determinant is 


’ ’ 2 - 2 A 2 
(23, - Zio)" = (MyyMo9 - Aygo, )"(Zaq > Zy9) (31) 
Substituting equations (29) and (31) into equation 


(25) shows that 


U(Z') = U(nZny) = U(Z) (32) 


20 





5. Example of a transformation to a unilateral form. 

The non-unilateral amplifier of Figure 8(a) can be 
transformed into a unilateral form rather simply. The 
voltage represented by I5Rj5 can be subtracted from the 
input terminals by means of a lossless transformer. Thus, 


Vz, in Figure 8(b) becomes, by proper choice of turns ratio, 


R 
Vi =v, - tev, (33) 
Roe 
so that the characteristic equationg of the transformed 
amplifier are 
vi - R41, + Ral Ri2(R I, + RooI5) 
E-~ ~li-l 12°22 © w—\"01-1 22-2 
R22 
R - Rank 
= %11R5 yoko, 
——S oo ey (34) 
R 
22 
V5 = (Roy ~ Rya)Iy + Reale 
or the transformed impedance matrix, Z', becomes 
RiiRoo ~- RyoRoq P 
ome Roo (35) 
Rai ~ Ryo Roo 


The requirements for stability under this transforn- 
ation are 

Riz eon 

Roo 7 O 
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Figure 8 


Rials 
= _ 


ee: 





(a) A non-unilateral amplifier 


ee Roo 
ToRi9 TyRoy 
+ * 
RyaiRoa 
€y! ——-°— y!-—> 
| ay Va 7” I, 


(b) Transformed Amplifier 


| ? 
Ry 2 


1 ap 
Fin t 
| — 
? q 
ty) b<—vj— 4—v3—» T2/ 
(c) Equivalent Unilateral Amplifier 
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6. Conclusions 

Every consideration of the stability of a network 
has as its basis the fact that the roots of the circuit 
determinant must lie in the left half of the complex fre-~ 
quency plane, The investigation of stability, then, requires 
an investigation of the roots of A by one meane or another, 
whether it be by calculations which may often be tiresome 
or impossible, or by means of the Nyquist diagram which sig- 
nals the appearance of any roots in the right half-plane 
without their actual location being neceasary. 

The use of matrices in the first of these methods 
provides a powerful means of formalizing and organizing 
the investigation, even though the reduction of actual 
calculations may be minor. Honnell's criterion, for in- 
stance, signifies a great reduction in calculations once 
the overall] transmission matrix is attained. Obtaining 
the overall transmission matrix, however, may require con- 
Biderable calculation, so that the overall work may not 
be reduced with respect to any other method, The advan. 
tages of using matrix analysis then lies in the organizing 
of the problen. 

The criteria discussed previously all have the dis- 
advantage of the type of answer they give -- the circuit 
either is, or is not, stable. If a circuit lies on the 


verge of instability, there are generally a large number 
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of parameters which may be varied to secure a more stable 
network. The designer, then, usually has a choice in the 
number of factors he may vgry, but at the same time he may 
alao be faced with the problem of interaction between these 
parameters which might be the cause of some unwanted effect. 

Essentially one might say, then, that the designer is 
faced with another design problem with the previous one 
acting only as a general basis from which to start. 

The unilateral ampi4fier has been chosen in this paper 
as the basis for the new design problem because the con- 
ditions for absolute stability are not involved. Matrix 
algebra helps in transforming a non-unilateral amplifier 
into this unilateral form. Although its use does not ac- 
tually synthesize the desired clrouit, matrix algebra may 
act as a strong guide in determining what a transformation 
does, and lead our'intuition to the network which will ac- 


complish our aim. 
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APPENDIX I 


The impedance transformation does not always involve 
inversion of Z. Mason / 13_/ demonstrates that by means of 
another transformation of Z' into Z" of the same general 


form as equation (17) 


Zz" = jy" a n'(Z' + jy')ut ny 


n'(Z' = gx’ + jx’ + jy')74 n, 


n! [ n(z + jx)71 n, + qx! + jy'| 32 ny, 
Choosing x' + y' #0, 


Z" ” jy" 


n' | n(z + jx)ci n, | =1 ny 


n'nw+(Z + jx)n, n, 
which has the general form 


Z" = jy" 2 n(Z + Jx)n, 


~ 
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APPENDIX II 
FEEDBACK COUPLING 


Once an amplifier is reduceq to a normalized, unilater- 
al, and hence positive-definite form, an additional lossless 
coupling, S, may be considered for feedback, and analyzed. 

It 18 convenient to visualize the S terminals shown in 
Figure 9 as transmission lines of characteristic resistance 
unity. On icles sie the incident waves S and reflected 
waves W' are related by the scattering matrix S /10_/: 


Wi S11 Sjo S13 Sy4\ / Wi 
W2 : 321 “820 825 “Say | | Wo (28) 
W5 S31 S32 533 S34] | Ws 
Wi, Say- Sao Syz Say] \ Wa 


Since we postulate matched loads to the lines, waves 


W, and Wy must ne,essarily equal zero, and S becomes 


Si, 0 83 0 
5., 0  ge= 40 

‘ ; 21 23 (29) 
S31 0 S33 0 
Sy, 9 Syz 0 


From equation (28) we can therefore extract the submatrix 


eg 





t 
(") "ar “as\ [Ma (30) 


t e 

“Oy! a is 

SOURCE is 2 
<w! w! ae 

LOAD 1 B heel 
> @2 
Figure 9 - Feedback coupling 
or, more concisely, 
W' s SW (31) 


For the greatest efficiency of tranamission, we let 
the reflected waves W) and Wz vanish, The goupling circuit 
is lossless and waves W and W' must carry the same power, 


or, 
Jwole + |wyle = [WP + [ws] ° (32) 
In matrix form, this is written | 
se Wy 
(ws wh). re (Ww, Wa} - a (33) 
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or, more concisely, 
to wht _ # 
—) oe (34) 
where * denotes the complex conjugate. 


However, from equation (31) 


4 ¢% = he fi 
We ow'* > W,S,8"W (35) 


where Wi = W,S,, according to the reversal rule //11_/, 50 
that the power conservation condition of equation (34) 
requires 

S;s* = I, the appropriate unit matrix 
or 


Sy = en (36) 
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